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UTIAS TECHNICAL NOTE NO. 260 


ERRATA AND ADDENDA 


Please replace pages 8 and 9 with the attached revised versions. 

In the original page 9 the remarks concerning integration of dv r and 
dv^ were inadvertent and incorrect. Valid integrations apply only to the 
Cartesian forms du, dv, dw as set forth in Eq. (la) of the new page 8. The 
components dv r and dv^ in cylindrical pclars serve only to delineate the 
transfer functions of the shock-interaction process in simplest form. In 
the revision, the opportunity has been taken to expand page 8 by providing 
equations that, hopefully, clarify the verbal statements. 


Page 5. Definitions of: 

_K: In more conventional notation "0" would be replaced by " 9 + 

%/ 2": see remarks below re definition of 9. 

r, 4> , X^: Delete "of field point". [Refers to components of _K and 

d_Z describing shear wave (Figs. 1, 2).] 

X_: Delete "or , r, <S> , Fig. 2)". 

0: 6 is the polar angle of the velocity vectors, proportional to 

d_Z, in the shear wave, Fig. 1. The polar angle of the wavenumber 
vector _K, normal to this transverse wave, is 0 + it/2. 


May 21, 1986 
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Abstract 


This work is a small extension of our NACA studies of the early fifties 
that predicted amplification of turbulence on passing through a shock wave 
(observed for turbulent boundary layers), as well as the generation of 
intense noise (observed for supersonic jets). The first solved the basic 
gasdynamics problem of the interaction of an infinite planar shock with a 
single three-dimensional spectrum component of turbulence (an oblique 
sinusoidal "shear wave"). The second developed the comprehensive 3D 
spectrum analysis necessary to generalize the scenario to the interaction of 
a shock wave with convected homogeneous turbulence. Numerical calculations 
were carried out to yield curves (vs. Mach number) of rms sound pressure, 
temperature fluctuation, and two components of turbulent velocity downstream 
of the shock, for two cases of preshock turbulence. The present numerical 
study reproduces these for one case and provides in addition their 
one-dimensional power spectra (vs. wavenumber or frequency). Ratios of the 
several postshock spectra to the longitudinal preshock turbulence spectrum 
(ID) have been computed for a wide range of Mach numbers; curves vs. 
wavenumber are presented for two scenarios of preshock turbulence: isotropy 

and axisymmetry, both based on the von Karman 3D spectrum. 
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Nomencl ature 


a = 1.3390, pure number 

B = 55/18-rca, pure number 

c^, c = pre- and postshock sound speeds, respectively 

c* = critical sound speed 

= flux of acoustic energy from unit shock area 

IjUrb = flux of preshock turbulence kinetic energy into unit shock area 

K = nondimensional wavenumber vector (K 15 K 2 , K 3 or K, 0, <j>, Fig. 1; 

also K r , K x , 9, Fig. 2) 

k_ = K/aL, actual wavenumber vector 

L = longitudinal scale of turbulence 

= U^/c^, preshock Mach number 

M 1 = U/c, postshock Mach number 

P = transfer function relating dp" to du 

P = ambient pressure 

p" = pressure perturbation 

I = transfer function relating dV to du 

I = ambient temperature 

r, X x = cylindrical coordinates of field point 

Uft» U = pre- and postshock stream velocity, respectively 

u, v, w = nondimensional components of velocity perturbation in directions 
X ls X 2 > X 3 , respectively (actual components/c*) 

X_ = nondimensional field point vector (X x , X 2 , X 3 , Fig. 1, or X p r, 

4>, Fig. 2). X x is 1 shock 

x^ = J(aL, actual field point vector 

X = transfer function relating du' to du 

Y = transfer function relating dv' r to dv p 
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r 

9 


p A» p 
t' 




a 


a) 

Note: 


= (complex) wave amplitude associated with a= u,r, <t>, %' , p :i , 
etc. 

= ratio of specific heats (taken as 1.4) 

= polar angle (Figs. 1-3) 

= critical angle separating evanescent and nonevanescent pressure 
waves 

= fluid density upstream and downstream of shock, respectively 
= temperature perturbation 
= azimuth angle (Figs, 1, 2) 

= one-dimensional (ID) spectral density associated with 
a = u, v , u 1 , x‘ , p" , etc. 

= radian frequency (2 ^frequency ) 

superscript " signifies actual dimensional quantity. 
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INTRODUCTION 


Analytical studies of three-dimensional disturbances convected through 
— and interacting with — a shock wave seem to have commenced in the early 
1950's. 4 The disturbances were waves of (1) vorticity, (2) entropy, or 
(3) sound; they were in the form of oblique planar, usually sinusoidal, 
patterns. Any one such wave encountering the shock would generate all three 
kinds on the downstream side. It was pointed out by Ribne^’ 5 ’® and by 
Batchelor' that the vorticity waves (called therein "shear waves") were 
three-dimensional Fourier components of arbitrary incompressible flows, 
e.g., weak turbulence. (See also Moyal . ) 

The initial papers^ -4 treated the interaction of individual waves of 
arbitrary inclination with an infinite plane shock; the analyses were 
linearized in terms of wave amplitude. A later paper® developed the 
comprehensive 3D spectrum analysis necessary to describe the interaction of 
turbulence with a shock; the earlier single-wave results were the "building 
blocks". Numerical calculations were carried out for rms values of 
turbulence velocity components, temperature, and pressure (sound) 
fluctuations downstream of the shock, assuming either isotropic or 
axi symmetric preshock turbulence. The plotted curves showed that velocity 
components of the post-shock turbulence would be amplified as much as 45% 
relative to preshock levels. Recent measurements do display the phenomenon 
of amplification (e.g.. Refs. 9, 10). 

The calculation of such turbulence amplification on passaae through a 
shock has received renewed attention in recent years . 1 The term 
turbul ence in this context is, however, a misnomer: these papers have 
reverted to dealing with a single 3D spectral component of turbulence. They 
oversimplify in interpreting the single-wave results as representative of 
the broad spectrum of waves constituting turbulence. On the other hand, two 
of thenw^’ 1 ^ pioneer in the application of numerical integration of the 
Euler equations to the single-wave shock-interaction problem. By coping 
with nonlinearity they test the range of applicability of the results of the 
linear analysis. 

Some other studies based on the same fluid dynamics, but quite distinct 
in orientation, may be noted. These relate to the passage of a columnar 
vortex broadside through a planar shock wave: a cylindrical sound wave 
appears on the downstream side, partly cut off by the shock, as well as a 
modified vortex. The single-wave (single Fourier component) results of, 
e.g.. Ref. 1, can serve as the "building blocks" to compute this sound wave 
in detail. 4,15 Unlike the shock-turbulence interaction problem, which is 
stochastic, this phenomenon is deterministic (and is more readily verifiable 
by experi ment 1& » 1 ' ) . Nevertheless, the underlying analytical framework is 
the same. Both problems have been considered relevant to the generation of 
"shock noise" by turbulence passing through the shock pattern of a 
supersonic jet. 

It appears then that, despite other results both old and new, the only 
genuine calculations of turbulence interacting with a shock wave are those 
of the early reference 6. But these, as noted above, have been limited to 
rms values of the relevant disturbances: there are no computed spectra. 
The present study is an extension of that paper (in a very limited sense); 
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it seeks to provide the one-dimensional power spectra (vs. wave number or 
frequency) of velocity, temperature, and pressure perturbations, and of the 
acoustic energy flux emanating from the shock. The procedure is one of 
numerical integration of the corresponding 3D spectra. These are adapted 
from Ref. 6, with the 3D spectrum of the preshock turbulence specified. 


SHOCK-TURBULENCE INTERACTION 


Transfer Functions Connecting Fourier Components (Deterministic) 

A snapshot of an arbitrary flow field may be represented by a Fourier- 
Stieltjes integral in three dimensions, as 

u(X) = Jdu = JdZ u (K)exp(iK-X); K = K x , K 2 , K 3 

v(X) = Jdv = JdZ v (K)exp(iK/X_) ; |J<| = K = 2-n:/wavel ength (la) 

w(X) = Jdw = /dZ w (K)exp(iK-X); X = X x , X 2 , X 3 

This integral is effectively a superposition of plane sinusoidal waves with 
wavenumber J< normal to the planes of constant phase; variation of J< implies 
a distribution of wavelengths and ori entati ons . We can apply this format to 
weak, essentially "frozen" turbulence [a pattern with negligible time 
dependence); this will behave almost incompressi bly , governed by 

su/aX! + av/ax 2 + aw/ax 3 (lb) 

even though convected at high speed. Applying this constraint to (la) shows 
that the amplitude dZ and wave number K_ are orthogonal: 

K x dZ u + K 2 dZ v + K 3 dZ w = K°dZ = 0; dZ_ = dZ u , dZ v , dZ w (1c) 

Thus an individual wave is transverse; it may be interpreted physically as 
an oblique sinusoidal wave of shearing motion (Fig. 1). 

Such a wave, when convected into a shock, interacts in a predictable 
fashion according to linear theory. ^ a "refracted" shear wave, a 
superposed entropy wave, and a pressure wave emerge on the downstream side. 
If the initial pattern of waves (upstream turbulence) is known only 
statistically, then the downstream pattern (modified turbulence, entropy 
"spotti ness" , and noise) can be determined statistically. That is, spectra, 
correlations, and mean square values can be calculated. 

1 

To this end, we incorporated the deterministic single-wave relations 
into a comprehensive spectrum analysis in three dimensions” for homogeneous 
turbulence. A brief account of relevant parts of the development is given 
below. The physical quantities are normalized so as to be nondimensi onal : 

u, v, w, etc. = velocity components/critical sound speed, c* 
p" = pressure perturbati on/ambient, P A 
x* = temperature perturbation/ambient, T 

(But addition of a superscript ~ to u, p" , %' , etc., removes the normaliza- 
tion.) 



It will be convenient to re-express the velocity field of an incident 
shear wave (Fig. 1) in cylindrical coordinates as (see also Fig. 2 where the 
wave is viewed edge-on) 


du = dZ u exp(i K/X) 

dv r = dZ r exp(i K_*X) (2) 

dv* = dZ, exp(i K*X_) 


where du is parallel to Xj (normal to the shock), dv r is parallel to r, and 
dv^ is perpendicular to r and in the direction of increasing <t>. The 
planes of constant phase _K*_X (= _k*x, see below) = constant make an angle 0 
with the X x axis. The wave number K and position vector X are non- 
dimensional; they are formed from their dimensional counterparts ]< and x as 


J< = j< a L ; L = turbulence longitudinal scale 

X, = x/aL; a = pure number 

respecti vely . 


(3) 


Figure 3 shows the results of the encounter of the incident wave, Eq. 
(2), with the shock. The three waves that appear on the downstream side 

are: 


refracted shear wave , with components: 

du' = dZ^(J<' )exp(i _K' vX) ; 
dVp = dZp(J<' )exp(i JK'*X.); 
dv 4> = dZ^JC' )exp(i K' *X); 


dZ u = * dZ u 
dZp = YdZp 

= dZ^ 


entropy wave , aligned with refracted shear wave: 

dT 1 = dZ! c (K , )exp(i K' -X); dZ^ = TdZ u 


sound wave : 

dp" = dZp(J<" )exp(i K" *X); dZ|J = PdZ u 


(4) 

(5) 

( 6 ) 


The planes of constant phase _K* »_X = constant and K" *X = constant make angles 
0' and 0", respecti vely, with the X^ axis or shock normal (Fig. 3). These 
angles are functions of the angle 0 of the incident wave (and of the Mach 
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number). For |0| in the range from 'ZercT'to a critical value 0 cr , the 
pressure wave decays exponentially with distance from the shock (in 
proportion to wave number): such waves are called evanescent . For 1 0 1 in 
the range from 0 cr to n/2 the wave shows no decay and is called 
nonevanescent herein. 

The right hand set of equations relates the respective amplitudes of 
the downstream waves with those of the initial shear wave components, dZ u , 

dZ p , dZ.. The transfer functions X, Y, T, and P are all dependent on the 
incident wave angle 9, as well as Mach number M; they are the results of a 
gasdynamic analysis carried out in Ref. 1. Functional expressions, taken or 
adapted therefrom, are given in Appendix A herein. The transfer functions 
are quite different in the regimes 0 to 9 cr and 0 cr to n/2; in particular, 
the form of P dictates evanescence in the former and nonevanescence in the 
latter . 


Three-Dimensional and One-Dimensional Power Spectra (Stochastic) 

For application to a stochastic field such as turbulence it is 
necessary to go over to statistical relations. If we form an ensemble 
average for waves of different wave numbers _K and £ in respective ranges d 3 K 
and d 3 < it is easily proved that 


dZ u (K)dZ*( £ ) = [uu] &(< - K)d 3 K d 3 £ , (7) 

if the turbulence is homogeneous. 8519 The quantity [uu] is a special symbol 

for the three-dimensional spectral density of u 2 in wave number space _K. 
The vanishing of 6(jc-K) for £ f K_ implies that waves of different 
wavelengths or inclinations (since £ and _K are vectors) are statistically 
independent. 

The integral of (7) over £-space may be written, by virtue of the 
5-function, 


dZ u (K)dZ*(K) = [uu] d 3 K (8) 


In similar fashion we can develop corresponding expressions for the 3D 
spectra [u'u 1 ], [t't 1 ], [p"p"], etc. The integral of (8) over _K-space is 

u 2 , that of [tV] over K'-space is x' 2 , and so on. 

Application of these procedures to (4) to (6) leads to 6 
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(9) 


u 2 = /// [uu]d 3 j< 


U ' 2 = III [u‘u*]d 3K' 


T* 2 = III [x'x'ld^K 1 


p" 2 = ffj [p"p"]d 3 j<" 


» /// I X| 2 Cuu 3 d 3 K_ 
= Ilf | T | 2 [u u ]d 3 j( 

= HI | P| 2[uu]d 3K 


where the integration limits are ± °° in K t , X 2 , K 3 , etc. 

We shall limit attention to axisymmetric turbulence, with main emphasis 
on the special case of isotropy. This and the axisymmetry of the shock 
interaction process led us in the earlier material” to employ cylindrical 
coordinates. In these coordinates. 


d 3 K = d<|> K r dK p dK 1 (10) 

A first integration with respect to d<t> then yields a factor 2n, so that 


u 


2 


u 


» 2 


x 


• 2 




2tc I 
2n j 
2n j 


I" [uu]K dK dK x = /“ 
o -co 


r |X|2[uu]K dK dK x = 
o 


f | T| 2 [uu]K r dK r dK 1 = 
0 


f I P| 2 [uu]K r dK dK x = 
o 




J (K l )dK l 

— CO 


I ^ x .(K 1 )dK 1 

— CO 


r ® p "(K 1 )dK 1 

— CD ' 


( 11 ) 


respectively. 

The last equality in each line defines a one-dimensional spectrum; 
these are explicitly: 
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00 


4> u (K x ) = 2n f 

[uu]K p dK r 

0 


cf 
*— * 

II 

rv> 

' ‘ 8 

| Xj 2[ uu ]K r dK 

0 


CO 

(K x ) = 2n f 

1 T| 2[ uu ]K r dK 

0 


CO 

§pM (KJ = 2% l 

|P| 2 [uu]K r dK, 


0 


As formulated, all four spectra are functions of the upstream longitudinal 
wave number K x . This is convenient and puts them on a common basis. 
Moreover, K x is proportional to radian frequency to in the form 

K i = waL/U A ( 13 ) 

[This follows from 

kj = to /U A (14) 

and the definition, Eq. (3), connecting the two k's.] 

Although th§ interaction process of Fig. 3 conserves K r but not K,, 
(K x / K x f K x ) , the frequency to is invariant (to = to 1 = to"). 
That is, an observer perceives the same frequency as each of the wave 
patterns moves past. This is, of course, a necessary result for a 

statistically steady process. A geometric analysis formally confirms the 
invariance of to. 


Isotropic Preshock Turbulence 


For evaluation of the ID spectra of (11) the 3D spectrum, [uu], of the 
input turbulence must be specified. The von Karman spectral model (called 
9 11 in Ref. 20) is chosen; in our notation it is 


[uu] = u 2 


BK r 2 

2tc[ 1 + K l 2 + K r 2 ]17/6 


(15) 


B = 55/18*?, a = 1.3390 


where the longitudinal scale L of the turbulence is incorporated in K = kaL. 
The pure number a is a normalizing constant chosen so that 
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/ dK, / [uu]27iK r dK r = u 2 

-CO 0 


For numerical evaluation of the ID spectra it will be convenient to 
express the cylindrical wave number component K r in terms of the 
longitudinal component Kj and the polar angle 0. This will replace the 
infinite range in K r by a finite range in 0. It will also be easier to 
interpret the integral in terms of the geometry of Fig. 3. Compatibly with 
Ref. 6 (Eq. 56) we put (Figs. 1, 2): 


K 1 = -K si n 0, K 2 = K cos 0 cos <J> 


K 3 = K cos 0 sin*, K r = /K 2 2 + K 3 2 = jKjcote (17) 

dK r = \K X \ esc 2 0d 0 

Inserting (15) and (17) into (11) gives* after some reduction, the following 
format for evaluation of the ID spectra: 


where 


®i^l) = B f 71 ^ I r I 2 COS 3 0d0 

-^T“"k7^ J o 111 iil^0[b ir TT^0F 7 7^ 


(18) 


b 1 = (1 + K 1 2 )/K 1 2 


and 


for i = u, u 1 , t 1 , p" , r 3 


X T, P. 


These express the general form to be numerically integrated for the shock- 
interaction products in the present case; namely, isotropic preshock 
turbulence with statistics described by the von Karman 3D spectrum. 


One Dimensional Spectra of v 2 and v‘ 2 


The one-dimensional spectra of v 2 and v' 2 involve a less straight- 
forward derivation. The following section sets forth the basic equations 
and final spectra. 

For isotropic and axisymmetric turbulence, v 2 = w 2 upstream of the 

shock and v' 2 = w' 2 downstream of the shock: the respective ID spectra have 

similar equalities. The corresponding 3D spectra, [vv] and [ww] , [v'v‘] and 
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[w'w 1 ], are not respectively equal, nor are they axisymmetri c. The sums, 
however, can be shown to be axisymmetric : they depend on K r rather than 

on K 2 and K 3 separately. These sums on both sides of the shock have a 
relatively simple connection. 


[v'v 1 ] + [w ! w‘] = (I?! 2 - 1) tan 2 e [uu] + [vv] + [ww], (19) 


where the analytical form of the transfer function Y is given in Appendix A 
herein. 

We now restrict ourselves to the special case of isotropic preshock 
turbulence. Then invoking the cited equalities and the axisymmetry of (19) 
and employing the input spectrum and methods of the last section leads to 
the desired lD-spectra (details are given in Appendix B). The results are: 


V K l) _ B J 71 / 2 (2 + cot 2 9)cos9d8 _ W 

u 2 2K^ /3 0 si n 3 0[b' + cot 2 e] 17/6 u" 2 


( 20 ) 


V( K i) _ VKl) 


2K 


B 

"57T 


/ 


% ! 2 ( 1 Y [ 2 - l)cos9d9 

s i n 3 0[b ' + cot 2 0] 17/6 


V( K l) 


( 21 ) 


One-Dimensional Spectrum of Acoustic Energy Flux 


The last of equations (18) gives the ID spectrum of the sound pressure 
generated downstream by the passage of isotropic turbulence through a shock. 
The connection to acoustic energy flux is not trivial, as in the case of 
quasi -plane waves propagating through still air. An analysis has been given 
by Ri bne r^ j; i nv ok i ng energy flow relations for a moving medium 
( B1 okhi ntsev V . The integral of Eq. (15) of Ref. 21 reads, in the present 
not at i on , 

I AC = J— / dp" 2 ( 0“ ) (1 + M 1 sin0")(M 1 + sin0") (22) 

pc e cr 


This evaluates the flux of acoustic energy emanating from unit area of the 
shock. Since only nonev anescent waves figure in I A q, the limits of 
integration correspond to that range. The result after numerical 
integration is plotted in Ref. 21. 


The quantity dp" 2 (0) can be identified with the integrand (including 
the 2n factor) of the last of Eqs. (11) herein. Hence the further 
development leading to the last of Eqs. (18) may be applied. This yields 
the one-dimensional spectrum of this flux, ®AC ( K 3. ) , i n the form 
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9 


£ac (K i) = 2 j 2 B i P | 2 ( 1 + H 1 sin0")(M 1 + sin0" )cos 3 0d 9 

*TURB '^ 2 C * e cr sin 5 0 [b ' + cot 2 0] 17/6 


when the input turbulence has the von Karman (3D) spectrum. This flux has 
been ratioed to the flux of kinetic energy of turbulence through the shock. 


*TURB = f PA U A C * 2 u2 = | P Uc * 2 u2 ( 24 ) 

(The — perhaps surprising — factor 5/2 results from the definition 

Z TURB = j PA t( U A + c * u ) 2 + ( c * v ) 2 + (c*w) 2 ](U a + c*u) - I p A U A 3 (25) 

and the postulated isotropy, u 2 = v 2 = w 2 .) 


One-Dimensional Spectral Ratios; Lack of Uniqueness (Aliasing) 

The component of turbulent velocity normal to the shock — called u 
herein — has a central role in the analysis. The three dimensional 

spectrum of u 2 governs the shock-interaction effects, and the resulting ID 
spectra are conveniently scaled to u 2 . Suppose, however, that we scale 

these ID spectra to the ID spectrum of u 2 : then this ratio in each case may 
be regarded as a sort of power spectrum transfer function; it will be a 
function of K 1 . Up to this point these quantities have been evaluated only 
for the von Karman 3D preshock spectrum, a particular case of isotropic 
turbulence. The question may now be asked, how much will these spectral 
ratios change with changes in the preshock spectrum? 

Such changes will be explored in the present section. They are 

expected to be nonzero because the ID spectrum of u 2 (in K^ does not 

uniquely determine its 3D spectrum (in K 15 K 2 , K 3 ) on which the shock 

interaction depends. (The 3D spectrum, on the other hand, does uniquely 

determine the ID spectrum as a result of the double integration over K 2 and 
l< 3 .) This lack of uniqueness has been termed an al iasing effect. We shall 
apply the term as well to lack of uniqueness in the postshock/preshock ID 
spectral ratios. 

A geometric interpretation may be developed with the aid of Fig. 2. 
The 3D spectrum of turbulence is an aggregate of waves like the one shown 
with a wide range of wavelengths 2n/K and angles 9 (and <t>) . If ®(K^ is the 
ID spectrum, the differential ®(K 1 )dK 1 is a narrow band of these waves with 
component wavelengths close to 2n/K 1 . It is evident that a variety of wave 
inclinations 0 with appropriate wavelengths 2-rc/K (functions of 0) could be 
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chosen to have this same x 1 -di rection intercept 2n/K 1 . Thus the 3D wave 
pattern corresponding to $(K X ) is not unique. 


Axi symmetric Turbulence 


A convenient example of aliasing changes in spectral ratios can be 
obtained by modifying the preshock spectrum to change it from full isotropy 

(v 2 = w 2 = u 2 ) to axial symmetry (v 2 = w 2 =/ u 2 ). The conditions for a 3D 
spectrum tensor to be axi symmetric have been discussed by Batchelor^ and 
others; for full generality they are very complex. However, we are 
concerned here with but an example of axi symmetry. This can be obtained 
from the von Karman spectrum (a particular case of isotropic turbulence) by 
a simple modification: we merely multiply the longitudinal spectrum [uu] 

(the only one of concern herein) by an arbitrary non-negative function of 
wave angle 0 (or of K]/K r = tan0). (9 is considered to be a polar angle 
and is restricted to the range 0 to ix/2.) 

This axi symmetric spectrum then takes the form 


[uu] AX i = F 2 ( 0) [uu] 


(26) 


This is justified by comparison with the second of Eqs. (12) 


[uV] = | X| 2 [uu] 


(27) 


that relates the postshock longitudinal spectrum, [u'u 1 ], to the preshock 
value, [uu], X being the transfer function. There being no preferred 
direction in the plane of the shock (1 to u) , [u'u‘] must surely be 
axi symmetric. But the argument for axi symmetry does not depend on the 
particular form of the factor X: it could be an arbitrary function of 

Kj/Kp (or of 0). 

For our particular example we shall take 


[uu]axi = [uu]cos N 0, with N = 2 


(28) 


to describe the 3D axi symmetric preshock longitudinal spectrum in terms of 
an isotropic form [uu]; specifically, the latter is taken as the von Karman 
form that we have used throughout. For this case the ID spectra 
$ u i(K i ), $ 1 .i(K 1 ), and $_n(K 1 ) are given by Eqs. (12) with an extra factor 
cos 2 0 in the integrands. v 
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RESULTS AND DISCUSSION 


Isotropic Preshock Turbulence 


RMS Values of Postshock Disturbances , and Noise in Decibels 


Figure 4 gives the variation with upstream Mach number of the various 
shock interaction products for a specific scenario: the preshock turbulence 

is isotropic with an intensity of 1% of free stream. The curves display rms 
perturbations of longitudinal velocity u and lateral velocity v or w in 
percent of initial stream velocity, and of rms temperature and pressure 
(noise) in percent of ambient. The figure is adapted from Fig. 4 of Ref. 
6. The curves represent, in effect, the integrals with respect to K x of the 
respective one-dimensional spectra; that is, the integrals displayed in Eqs. 
(11). (The actual procedure, however, bypassed the ID spectra and employed 
only the specification of preshock isotropy. The results are independent of 
the preshock spectra, 3D or ID, so long as they are consistent with 

isotropy. ) 

The amplification of both the longitudinal and lateral components of 

the postshock turbulence is evident; it reaches some 45%, as noted in the 
Introduction, for the lateral component at high Mach number. The other two 
curves in Fig. 4 refer to the temperature and pressure fluctuations, 
respectively. (These are both spatial and temporal : rms values are the 

same from either point of view.) In first order , these are absent from the 

postulated upstream flows (extremely weak second order pressure and 
isentropic temperature fluctuations are associated with the specified 1% 
preshock turbulence). 

On an acoustical basis the pressure fluctuation (noise) generated by 
the shock-turbulence interaction is very intense. This is shown in Fig. 5, 
where the noise level in decibels corresponding to Fig. 4 is plotted vs. 
Mach number. (The definition is 


dB = 20 log 10 (/7VP REF ); P REp = 2x10-1° atm (29) 

when the postshock ambient pressure is taken as 1 atmosphere.) For 1% 
preshock turbulence the postshock noise level is predicted to exceed 140 dB 
for all upstream Mach numbers above 1.05. 

Figure 6 displays normalized ( nondimensi onal ) ID power spectra 
calculated from the equations presented herein; the scenario is isotropic 
turbulence, governed by the von Karman 3D spectrum, being cohvected by an 
M = 1.25 flow into a normal shock. The six spectra are: 


longitudinal comonent of preshock turbulence 
longitudinal component of postshock turbulence 


V K i)/^ 
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a T ,(K l) /f2 
[Ip.. ( K 1 )/P 2 ^ x =o 


[Sp"(Ki)/P 2 ] x= 

^AC^l)/ I TURB 


temperature fluctuation 

pressure fluctuation just downstream of shock (acoustic 
near field) 

pressure fluctuation far downstream of shock (acoustic 
far field) 

acoustic energy flux emanating from shock 


The first two of these are normalized by u 2 , which is the integral (from 
-°° to 00 in K x ) of ^.(Kj) (thus the integral of the first is unity). The 
next three are normalized by one or the other of (ambient temperature) 2 and 
(ambient pressure) 2 . Finally, is normalized by the flux of kinetic 
energy of preshock turbulence flowing into the shock. The superscript 
signifies a dimensional value, unsuperscri pted forms having been 
nondimensional i zed at the outset of the analysis. 

The pressure field (noise) decays from an extremely high value just 
downstream of the shock (X=0) to an asymptotic lower value — still very 
intense acoustically — far downstream (X=°°). Figure 6 shows a major 
difference in their spectral content: the X-0 near field is dominated by 

low frequencies, decaying asymptotically like K, beyond K^S. The X==° 
far field is very deficient in low frequencies; on a linear scale it has a 
bell-shaped spectrum, peaking near K 1 =l, but with the same asymptotic decay 

(the Kolmogorov K x law) beyond K-^3. 

One-Dimensional Spectra and Postshock/Preshock Spectral Ratios 

Figure 6 applies for M=1.25; a series of such figures could be 
presented for a wide range of Mach numbers. A much neater alternative, 

however, is to ratio each of these spectra at each value of K 1 to the %/u 2 
spectrum. This ratio, as has been mentioned earlier, could be regarded as a 
sort of transfer function connecting the ratioed pair of spectra. In this 
format the variation with Mach number can be discerned much more 
systemati cal ly. 

Figure 7 presents such spectral ratios: it relates the postshock to 

the preshock longitudinal component of turbulence (the u 2 divisors cancel) 
for a series of Mach numbers, M. For convenience the curves are normalized 
by factors Z(M) (tabulated on the figure) to force agreement with the M=1.25 
curve at K 1 =l; this makes the family of curves much more compact. It is 
seen that increasing Mach number enhances the low frequencies of the 
longitudinal component of the postshock turbulence. 

Figure 8 gives the corresponding spectral ratios relating postshock 
1 ateral component of turbulence to preshock longitudinal component. Here 
the behaviour is the reverse of that in Fig. 7: increasing Mach number for 
the most part decreases the low frequency content. 
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Figure 9 gives the spectral ratios relating postshock temperature 
fluctuation (arising from entropy "spottiness" generated at the shock) to 
preshock longitudinal component of turbulence. It is evident that the low 
frequencies in the temperature field are somewhat enhanced compared with 
those of the turbulence field. The variation is not, however, monotonic 
with Mach number: there is a foldover of the curves with increasing M. 

Figures 10 and 11 give the spectral ratios relating the near field and 
far field pressure fluctuations (noise), respectively, to the longitudinal 
component of the preshock turbulence. The two families of curves are 
clearly very different. Moving on. Fig. 12 displays the spectral ratios 
connecting the acoustic energy flux emitted by the shock on the downstream 
side to, again, the longitudinal component of preshock turbulence. One 
notes the marked qualitative similarity with Fig. 11. This is not 
surprising, since acoustic energy flux and far-field mean square pressure 
fluctuation are closely related. For a medium at rest the two are exactly 
proportional, whereas in the postshock flow they differ as described by 
Blokhintsev^ 2 for a moving medium. 


Axi symmetric Preshock Turbulence 


Postshock/Preshock Spectral Ratios; Aliasing Effect 

Calculations parallel to some of those pertaining to Figs. 7 to 12 
(isotropic turbulence, von Karman spectrum) have been carried out for a 
specified axi symmetric preshock turbulence. As described in the relevant 
earlier section, the two 3D preshock longitudinal spectra are related by Eq. 

(28). 


Figure 13 gives the ID spectral ratio of postshock to preshock 
longitudinal component of turbulence for the axi symmetric preshock 
turbulence specdified above, when [uu] is the von Karman spectrum (see, 

e.g.. Ref. 20). Also plotted is the ID spectrum $g/u 2 of the longitudinal 
component of the preshock turbulence. (Normalization of $ u is by the same 

u 2 as for Figs. 7-12, not by u^j.) Comparison with Fig. 6 (isotropic 

turbulence) shows that the two spectra of $ u /u 2 are vastly different. 
Nevertheless, the spectral ratios of Fig. 7 (isotropic turbulence) and Fig. 
13 (axisymmetric turbulence) are qualitatively rather similar, although 
quantitatively different. The spectral ratios are, in fact, very much less 
dissimilar than the spectra of the two kinds of turbulence. This property 
points up the utility of the ID spectral ratio in describing 
shock-turbulence interaction effects. 

The fact that the spectral ratios do change (although not radically) 
with marked changes in the 3D spectrum of the preshock turbulence is a 
consequence of the aliasing effect discussed earlier. Comparison of Figs. 7 
and 13 gives a measure, for a particular pair of cases, of this aliasing 
effect . 

Figure 14 gives the spectral ratio relating postshock temperature 
fluctuation to preshock longitudinal component of turbulence. This figure 
(axisymmetric preshock turbulence) is to be compared with Fig. 9 (isotropic 
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turbulence). Again, the sets of spectral ratios are qualitatively generally 
similar. The differences (aliasing effect) are greater than those between 
Figs. 7 and 13. 

Finally, Fig. 15 gives the spectral ratio relating the far field (X=») 
pressure fluctuations (noise) to the longitudinal component of preshock 
turbulence. This figure (axisymmetric preshock turbulence) is to be 

compared with Fig. 11 (isotropic preshock turbulence). The same sort of 
qualitative similarity is found between Figs. 9 and 14. The first glance 
impression of the similarity is, however, a bit misleading, e.g., for K<1 it 
is the M=1.01 curve of Fig. 11 that is very similar to the M=1 . 05 curve of 
Fig. 15. 

These three figures, taken together, give some idea of the changes in 
spectral ratio (the aliasing effect) due to changes in the 3D spectrum of 
the preshock turbulence. But this one comparative example of axisymmetric 
vs. isotropic turbulence hardly quantifies the effect. 
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APPENDIX A 


GENERAL RELATIONS AND TRANSFER FUNCTIONS 

The upstream Mach number M and the incident wave inclination 8 (Figs* 
1-3) are specified. These dictate a virtual Mach number R. A number of 
general relations are independent of the magnitude of R; the remainder, 
notably the transfer functions , take different functional forms, depending 
on whether R<1 or >1. The ratio of specific heats y is taken as 1.4. (For 
formulas in terms of y see Ref. 6.) 


General Relations 


M 

= 

speci fied 

m 

= 

6M 2 /(M 2 +5) 

M i 

= 

/(M 2 +5)/(7M 2 -l) 

0 

= 

speci fied 

0' 

= 

tan -1 (mtanO) 

R 

= 

Mj/cos 0' 

P 2 

= 

1-M 2 


= 

/|1-W 2 | 

P 

= 

tan“ 1 ( 1/ 

9 cr 

= 

cos -1 M 1 

0 cr 

= 

cot -1 (mcot0^. r ) 

9 cr 

= 

9 cr ' 71/2 

9" 

= 

, r M, 2 tan8' .. 

-tan _i [— — ] for W < 1 

P 2 


= 9 1 - p. for W > 1 


Transfer Functions 

These are defined in terms of functions of other functions. For 

numerical calculation the sequence should be in reverse order to the 
listing. The relative ease of programming belies an appearance of 
complexity. 
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w < l (o < | 0 | < e cr 

W > ] 

■ (0 cr < |9| < n/2 

p 

= 2 . 8 n M J 5 p 

( 2. 4m-. 4 )cos 0cos 9' 

P 

= No Change 

n 

_ (cose) /c2+d 2 0 “ K r x Po/^ 2 

n 

cosesinp. 

= c 


m 0 

mcos 8" 

T 

= /(atan0-l) 2 + (btane) 2 

T 

= N.C. , with b = 0 


*0.8(m-l) 2 /[2.4m-.4)/m] e 1 5t 



6 S 

= tan _1 (-B/A) 

6 s 

= 0 

6 t 

= tan _1 [b/(cot0-a)] 

5 T 

= 0 

s p 

. c8 w - dtanO 1 
= tan" 1 w 

dp w + ctan0' 

6 P 

= 0 

X 

1 5c 

= [Scos 9 '/cos 0Je ^ 

X 

= N.C. 

Y 

i 6c 

= [Ssin9'/sin0]e ^ 

Y 

= N.C. 

S 

= (/A 2 + B 2 /m)cos0 

S 

= N.C., with B = 0 

A 

= sece 1 + 2(m-l)cose‘ 

A 

= N.C. 


+ (a/m)(m-l) 2 sin0' 



B 

= (b/m) (m-1) 2 sin 0' 

B 

= 0 

c 

= D’a/m - F' 

c 

= N.C. 

d 

= D'b/m 

d 

= 0 

a/m 

= (CE + DF)/(C 2 + D 2 ) 

a/m 

= (C‘ + GF * )/ ( E * + GD’ ) 

b/m 

= (CF - DE)/(C 2 + D 2 ) 

b/m 

= 0 

C 

= (1/6 + 2m/3)tan 0' 

C* 

= m/3 - 2[1 + (m-1 )cos 2 0' ] 


-(l/2)[(m-l) 2 +(m-l)/1.2]sin29' 



D 

= D'p w /p 2 

D 

= Not used 

D 1 

= (m-l)[l + (m-l)cos 2 0'] 

D' 

= N.C. 

E 

= 2 - m/3 + 2(m-l)(p w 2 /p 2 )cos 2 0‘ 

E* 

= F ' (m-1 )/2 - ( l+2m/3)tan 0' 
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F 


F 


Not used 


= F'j3 w /p 2 

F* = (m-l)si n20' F ! = N.C. 

G = Not used G = -tan 9" 

Note: Errors have been found in one or the other of the two versions of 

Ref. 6 from which the above are taken or derived; these occur in several 
equations and in the list of symbols. Those that are applicable have been 
corrected. Some other errata are corrected in Appendix D herein. 

APPENDIX B 

ONE-DIMENSIONAL SPECTRA OF T 2 AND v 7 ? 

The starting point is Eq. (19): 

{[v ' v 1 ] + [w'w 1 ] }d 3 K_ = ( | Y | 2 - 1) [uu] tan 2 0 d 3 K_ 

+ {[vv ] + [ww]}d 3 J< 

Invoking the axisymmetry, we may integrate as in Eqs. (11) and (12): 

2 it / {vV] + [w'w']}k dK = 2n J (|Y| 2 -1) [uu] tan 2 0 K dK 
o ' o 

+ 2 n J {[vv] + [ww]}K r dK r (Bl) 

o 1 

The left-hand side is the sum of <^ / i(K 1 ) and §> i(K^); they are equal s even 
though their respective integrands are not. This implies 

2V(Kj) = 2V(k 1 ) = 2 it f(m 2 -l) [uu] tan 2 0K r dK r 

+ 2it / {[vv] + [ww] } K r dK„ (B2) 

o 

Applying the same argument to the right-hand side shows that the second 
integral may be interpreted as 2^,(K 1 ) = 2^(K 1 ). 


A. 3 



Isotropic Turbulence 


If we designate u, v, w as q l5 q 2 , q 3 , then the general 
for isotropic turbulence is:' 


[q.q.] = I£J<1 (K26 ii - K^.-) 
J 4itK 4 J J 


We take E(K)/47tK 4 to have the von Karman form 


E(K) = B u 2 1 
4nK 4 Zn [1+K 2 ] 1776 


Then, after a little manipulation, 


[uu] = liif ^ 1 ! 

211 [1+K L 2 +K r 2 ] 17/6 


[vv] + [ww] 


2K, 2 + K r 2 

2n [1+K X 2 +K r 2 ] 17/6 


where K p 2 = K 2 2 + K 3 2 . Insertion into (B2) then yields 


■MM - w sL tC 


(2K 1 2 +K r 2 )K r dK r 

[1+K 1 2 +K r 2 ] 17/6 


= V( K i) = V K i> 


+ iZ r (1 ^ 2 “ 1 ) K r 2 ( tan 2 e ) K r dK r 
2 0 [1+K 1 2+ K p 2 ] 17/6 


spectral form 
(B3) 

(B4) 

( B5) 

( B6) 

(B7) 

( B8) 


Invoking polar coordinates in the form of Eq. (17) yields: 



(B9) 


M K i ) = M K i ) = B -- u2 J TC/2 ■ ^ + -^ 2 - B -) c ° - s6 d9 

2K 1 5/3 0 sin 3 0[b'+cot 2 e] 17/6 


VW = V( K l) = M K i) 


B u 2 j 71 / 2 ( [ Y| 2 -l)cos0 d9 

2K x 5/3 0 sin 3 0[b l +cot 2 0] 1 - 7/ ' 6 


(BIO) 


where b‘ = (1+K 1 2 )/K, 2 , and Y is defined in Appendix A. Equations (B9) and 
(BIO) are the desired one dimensional integrals [cited as (20) and (21), 
respecti vely, in the main text]. 


APPENDIX C 
OBLIQUE SHOCKS 


The analysis and results, on a wavenumber basis, may be applied to 
oblique shocks by the usual procedure . The equivalent normal shock 
transformation 


M = M qCos 4* or U A = U 0 cosc|, (Cl) 

is made, where M Q is the upstream Mach number, and 4> is the oblique angle 
between the shock normal and the upstream flow direction. In Figs. 4 and 5, 
the designation "1% preshock turbulence" is now to be interpreted as 
"(cosc|>)% preshock turbulence". 

The interpretation of all figures in terms of wavenumber l< x is 
unaltered, except for Fig. 12 (acoustic energy flux). That figure is 
inapplicable to oblique shocks: the structure of the underlying equations 

(22) and (23) is changed, and no simple proportionality will serve. 

The proportional ity of K 1 to frequency , however, does not carry over to 

the oblique shock case:* equations (13) and (14) are inapplicable. If K l 
is the component wavenumber along the stream velocity U Q , then the two 
equations are replaced by 


*An exception is the case of d? u ( K L ) for isotropic preshock turbulence. 



and 


= o)aL/U 0 

( C 2 ) 

= w/U 0 

( C3) 


respectively. Thus K l takes over the role of (and U 0 the role of U^). 

It follows that ID spectra in terms of K, , rather than K,, are 
equivalent to frequency spectra in the oblique shock case. This implies a 
coordinate axis rotation of amount 4 > about the 0X 3 axis: 


with inverse 


K x = i^coscj) + K 2 sin4) 

K 2 = -K x si n cp + K 2 coS( 1 j 

k 3 = k 3 

K l = K]COSc|> - K 2 sin4> 
K 2 = K x s i n cp + K 2 cos 4 j 


(C5) 


(C 6 ) 


Also spherical polar and cylindrical coordinates are defined about the 
K 1 -axis, similar to (17); thus (since the magnitude K = K) , 


= -Ksine, 


K 2 = Kcose cos 4 . 


K 3 = Kcos 9 sin<)> 



|K l |cot 0 


dK p = |K 1 |csc 2 9de 


( C7 ) 


In order to obtain the various ID spectra in terms of rather than 
K 1 , Eq. (12) must be rewritten with K p and (implicitly) <j>, 0 replacing K r , 
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9 S ^respectively. Also the factor 2% must be replaced by an integration 
over d4>; this is required since the transfer functions | r-j | depend on both 4 > 

and e, not 0 alone. Stated otherwise, the | r i | are axisymmetric relative to 
the shock normal direction (along which K x is measured), but not 

axisymmetric relative to the oblique-shock flow direction (along which K, is 
measured). The ID spectra are thus 


(K x ) 


/ d 4> / 
o o 


r.- 1 2 [uu]K r dK r 


( C8) 


Because of the postulated isotropy of the preshock turbulence, [uu] has the 
same form in the rotated reference frame as in the original frame. Thus it 

/v (V 

takes the form (15), but with K, , K_ replacing K ls K respectively. Then, 
with the use of Eqs. (C6) and (C7), the explicit form of the ID spectra is 


2n 


%/ 2 


= B_ r dp (• | r 12 

— ~5/3 J o_ J 1 

,, 2 K 0 0 


cos 3 9d 0 


si n 5 0[b 1 + cot 2 e] 17/6 


(C9) 


where \v^ \ = 1^1(0) is defined in Eq. (18), b' = (1 + K 1 2 )/K l 2 , and the 
angle 0 is evaluated from 


9 = tan- 1 (|K 1 |//K 2 2 + K 3 2 ) 


tan -1 


| cos - (K 2 /K 1 )sin<p| 


with 


[sin 2 cp + (K 2 /K j_) 2 cos 2 c|, + (K 2 /K 1 )sin2(J) + (K 3 /K 1 ) 2 ] 

K 2 /K x = -cote cos<l) 

K 3 /K 1 = -cote sin<j> 


172 


(CIO) 


(Cll) 


Equation (C9) , with (CIO) and (Cll), evaluates the ID spectra in the 
form of^a double integral. Thus two numerical quadratures are required (in 

<t> and 0) for the oblique shock case. By contrast, only a single numerical 
quadrature was required (in 0) for the normal shock case. 
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APPENDIX D 


CORRECTIONS TO THE BASIC REFERENCE, RIBNER ( 1954) 6 
The corrected equations are: 


— - — - | S | 2 sin 2 0‘ - sin 2 9 — — 

v 1 2 + w' 2 = J [uu] dk_ + v 2 + w 2 

COS 2 0 


[uu] = k -2 F(k)cos 2 0 


— go 2 n t e/2 

u 2 = 2 J F(k)dk / d<|> / cos 3 0de 
o o o 


2 % %/2 


u' 2 = 2 / F(k)dk / d <J> / |S| 2 cos 2 0'cos0d0 

o o o 


(54) 

(55) 

(57) 

(58) 


Corrections to the definitions in Table A (Symbols) are: 


Multiply by -1 


e" 


. M, 2 tan0‘ 

-tan -1 — , 0 < 0 < 0 r „ 


For e i n -re / 2 rea< j e -imn;/2 


There are a number of other typos that are either trivial or self-evident. 
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FIG. 1 PERSPECTIVE VIEW OF SHEAR WAVE (3D SPECTRAL COMPONENT OF 

TURBULENCE) IN RELATION TO REFERENCE FRAME. NOTE THAT d_Z IS NOT IN 
GENERAL IN THE PLANE OF K AND U. 



FIG. 2 PROJECTIVE VIEW OF SHEAR WAVE IN RELATION TO REFERENCE FRAME. 




FIG. 3 INTERACTION OF SHEAR WAVE WITH SHOCK: VIEW IN X x , 
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r-PLANE. 
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FIG. 4 AMPLIFIED TURBULENCE AND OTHER DISTURBANCES PRODUCED DOWNSTREAM OF 
SHOCK BY INTERACTION WITH ISOTROPIC TURBULENCE. PRESHOCK 
TURBULENCE INTENSITY 1 %. RMS VELOCITIES IN PERCENT OF INITIAL 
STREAM VELOCITY; RMS TEMPERATURE AND PRESSURE (NOISE) IN PERCENT OF 
AMBIENT. 
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FIG. 5 NOISE GENERATED BY SHOCK-TURBULENCE INTERACTION (ISOTROPIC 
TURBULENCE) ON A DECIBEL SCALE. POSTSHOCK AMBIENT PRESSURE, 
1 ATM. 





FIG. 6 NORMALIZED ID POWER SPECTRA OF THE SHOCK-INTERACTION PRODUCTS WHOSE 
RMS VALUES ARE GIVEN IN FIG. 4 (EXCEPT . /u*) PLUS TWO OTHERS (SEE 
TEXT). MACH NUMBER M = 1.25 ONLY. 
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FIG. 7 RATIO OF ID SPECTRA FOR ISOTROPIC PRESHOCK TURBULENCE: POSTSHOCK 

LONGITUDINAL COMPONENT/PRESHOCK LONGITUDINAL COMPONENT. 
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FIG. 8 RATIO OF ID SPECTRA FOR ISOTROPIC PRESHOCK TURBULENCE: POSTSHOCK 

LATERAL COMPONENT/PRESHOCK LONGITUDINAL COMPONENT. 
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FIG. 9 RATIO OF ID SPECTRA FOR ISOTROPIC PRESHOCK TURBULENCE: POSTSHOCK 

TEMPERATURE FLUCTUATION/PRESHOCK LONGITUDINAL COMPONENT OF 
TURBULENCE (BOTH NORMALIZED). 
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FIG. 10 RATIO OF ID SPECTRA FOR ISOTROPIC PRESHOCK TURBULENCE: 

"NEAR FIELD" NOISE (X = 0)/PRESH0CK LONGITUDINAL COMPONENT OF 
TURBULENCE (BOTH NORMALIZED). 
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FIG. 11 RATIO OF ID SPECTRA FOR ISOTROPIC PRESHOCK TURBULENCE: POSTSHOCK 

"FAR FIELD" NOISE (X = »)/PRESHOCK LONGTUDINAL COMPONENT OF 
TURBULENCE (BOTH NORMALIZED). 






FIG. 12 RATIO OF ID SPECTRA FOR ISOTROPIC PRESHOCK TURBULENCE: POSTSHOCK 

ACOUSTIC ENERGY FLUX (NORMALIZED BY PRESHOCK TURBULENCE ENERGY 
FLUX) /PRESHOCK LONGITUDINAL COMPONENT OF TURBULENCE (NORMALIZED BY 
MEAN SUARE VALUE). 
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FIG. 13 RATIO OF ID SPECTRA FOR AXI SYMMETRIC PRESHOCK TURBULENCE: 

POSTSHOCK LONGITUDINAL COMPONENT/PRESHOCK LONGITUDINAL COMPONENT. 
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FIG. 14 RATIO OF ID SPECTRA FOR AXI SYMMETRIC PRESHOCK TURBULENCE: 

POSTSHOCK TEMPERATURE FLUCTUATION/PRESHOCK LONGITUDINAL COMPONENT 
OF TURBULENCE (BOTH NORMALIZED). 
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FIG. 15 RATIO OF ID SPECTRA FOR AXISYMMETRIC PRESHOCK TURBULENCE: 
POSTSHOCK "FAR FIELD" NOISE (X = ~)/PRESHOCK LONGITUDINAL 
COMPONENT OF TURBULENCE (BOTH NORMALIZED). 





